[14] Particle Data Group, Phys. Rev. D50 (1994) 1472. 

[15] For details regarding renormalization on the lattice and lattice perturbation theory, see 
G.P. Lepage and P.B. Mackenzie, Phys. Rev. D48 (1993) 2250. 

[16] Although the value is not explicitly quoted, the world average which seems to have 
been used in the calculation is from the Particle Data Group’s 1986 number B 0 id(T —> 
v T 7r7r7r) = 8.3(8)% Particle Data Group, Phys. Lett. 170B (1986) 1. 

[17] One can update the phenomenological prediction using a simple scaling argument. Since 
the branching ratio scales with the square of the decay constant, we can shift the decay 
constant to what Isgur, et al., would have calculated with the present branching fraction: 

rExp _ rExp 

JCorrected IMR JlMR 

For example, a new world average of B new (T —> = 8.56(24)% [18] implies 

fCorrected imr = 0.254(20) which still agrees with the prediction of the a\ decay constant 
from lattice QCD. Note that even a 20% increase in the branching fraction will result 
in a 10% increase in the decay constant. This will only improve agreement between the 
experimental and theoretical decay constants. 

A similar statement may be made about a change in the world average of the r lifetime. 

[18] B. Heltsley, Tau 94/Montreux Proceedings. 

[19] This problem has a long and extensive history. Two recent papers are: H. Lipkin, Phys. 
Lett. B303 (1993) 119; M. Suzuki, Phys. Rev. D47 (1993) 1252. 

[20] For a discussion of these methods, see A. Duncan, E. Eichten, J. Flynn, B. Hill, G. 
Hockney, and H. Thacker, FERMILAB preprint FERMILAB-PUB-94.164-T; T. Draper 
and C. McNiele, Proceedings of Lattice ’93, p.453. 

[21] C. Bernard, J. Labrenz, and A. Soni, “Proceedings of Lattice ’92”, Amsterdam, 1992, 
Nucl. Phys. (Proc. Suppl.) B30 (1993) 465; Phys. Rev. D49 (1994) 2536. 
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/£•"' = fkr Mfe,)’ (i) 2 \ \ji “fr (1-0.31a,,), (14) 

where k c is the hopping parameter where the renormalized quark masses vanish; i.e. k —>■ n c 
is the chiral limit. 

With these caveats, we report lattice measurements of = 1390 ± 80 MeV and 

f Kl = 0.33 ±0.03 (GeV) 2 . 

We conclude by remarking that these calculations are extremely straightforward exten¬ 
sions of commonly performed lattice measurements. We encourage other lattice groups to do 
them, so that the sensitivity of these predictions to the standard lattice systematics (lattice 
spacing, simulation volume, quenching vs. dynamical quark mass) can be explored. We 
also encourage the lattice community to consider tau decay as a source of phenomenological 
problems. 
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that ay = 0.29(5) and 0.31(5). Therefore, the continuum a\ decay constant is calculated to 
be 


/« 


C ont 

ai 


0.33 ± 0.04 (GeV) 2 for am q = 0.010 


( 11 ) 


y 0.27 ± 0.05 (GeV) 2 for am q = 0.025. 

Again, we treat these results as independent determinations of the same quantity and average 
them to get 


fC m , _ o 30 ± 0 03 (GeV) 2 . 


( 12 ) 


The calculation of a branching fraction is very complicated due to the large width of the 
a\ [10,2], The value of f ai quoted in the 1989 paper by Isgur, Morningstar, and Reader was 
derived phenomenologically using the absolute rate for r —> v T tttttt [16]. They found that 
the experimental data at that time give the following decay constant [17] 

/Fmr = 0.25 ±0.02 (GeV) 2 . (13) 


This is in good agreement with our theoretical calculation. 

We can attempt a similar calculation of the mass and decay constant of the strange axial 
vector meson, the Ah- However, this calculation is not complete. The physical Jlh(1270) and 
Ah (1400) mesons are combinations of the 3 Pi and 1 Pi states [19]. While the correlator in 
equations (2) and (3) couples asymptotically to the lightest state which has nonzero overlap 
with its operators, the contamination of the next state only dies away as exp(—AMt), where 
AM is the mass gap. The two Ah states are close in mass, and since we do not have multiple 
operators to do a multistate fit [20], we cannot be sure that we are seeing the lighter Ah- 
Nevertheless, in the absence of other theoretical predictions, we attempt a calculation. 

Fixing the </ meson to its physical mass gives us /c s , the Wilson hopping parameter 
corresponding to the the strange quark mass. Then we are able to find the lattice mass 
and lattice decay constant by extrapolating the lighter quark to zero mass while holding 
the heavy quark at the strange quark mass. Finally, we must correct the factor of 1/4 in 
eqn. (10) since only one quark is light [15,21]: 
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(7) 


1270 ± 80 MeV for am q = 0.010 

M ai = 

1230 ± 140 MeV for am q = 0.025. 

Although the quantities we calculate on the lattice change with the sea quark lattice 
mass as seen in (5) and (6), the physical observables are independent of the mass of the 
dynamical fermions within errors. This has been true of all lattice calculations of the QCD 
spectrum to date. Therefore we treat the two simulations as independent and average the 
two results giving a physical a\ mass of 

M ai = 1250 ± 80 MeV. (8) 

This is in very good agreement with the experimental value 1230(40) MeV [14]. With these 
definitions of the lattice spacing the same simulations predict rho masses of 680(40) and 
640(60) MeV and nucleon masses of 990(60) and 920(80) MeV for lattice quark masses 
am q = 0.010 and 0.025, respectively. 

To compute the lattice decay constant, we perform a three-parameter correlated fit to the 
two different correlation functions for the three lightest quarkonium states and extrapolate 
to zero quark mass. We find 

' 

r 0.881 ± 0.036 for am q = 0.010 

/a = (9) 

0.803 ± 0.056 for am q = 0.025. 

In order to extract the continuum value from the lattice value, we must set the scale 
with the lattice spacing and multiply by a renormalization factor to convert from lattice to 
continuum regularization. We use the tadpole-improved formalism of Ref. [15], where 

fifi = /.“(aJf.j’Q 1 jd-0.31a,,). (10) 

and ay is the running coupling constant defined through the potential on the lattice at a 
momentum scale q = 1.03/a. Through the expectation value of the plaquette, we measured 
the strong coupling at q = 3.41/a and ran it down to q = 1.03/a using the two-loop QCD 
beta function. For our two values of the sea quark mass, am q = 0.010 and 0.025, we found 
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a\ in terms of the on-shell decay constant. The calculation is in complete analogy with the 
decay chain e + e“ — > p —> 7T7t, where the rho decay constant is taken to be a constant f p and 
the decay into pions is modeled by an on-shell rho decay matrix element. This calculation 
neglects processes where the a\ flows backwards in time. In addition to extracting f ai from 
data, IMR also present a calculation of it in a quark model; our work replaces the quark 
model calculation with a direct determination of the decay constant from QCD. 

Presently it is not possible in lattice QCD to include the decay of the a\ into p and 7r 
mesons, which would give us the off-shell form factors. 

In order to find the continuum a\ mass M ai and decay constant /^ ont , it is necessary to 
compute four quantities: the lattice mass aM ai , the lattice spacing a, the strong coupling 
constant defined in a particular prescription ay, and the lattice decay constant /^ at . 

The lattice mass is measured by performing a two-parameter correlated fit [11] to the 
axial-vector propagators for the three lightest quark-antiquark pairs. The results of these 
fits are presented elsewhere [12]. Since the valence quarks in the simulation are much heavier 
than the physical u and d quarks, we extrapolate linearly to zero quark mass. We find 


aM ai 


' 

0.670 ± 0.018 for am q = 0.010 
0.738 ± 0.049 for am q = 0.025. 


(5) 


The errors in these numbers are purely statistical. 

One expects [13] the mass difference between heavy QQ L = 1 states and L = 0 states 
to be weakly dependent upon the quark mass, and so we use the quarkonium S-P mass 
splitting to determine the lattice spacing. This calculation is performed in our companion 
work [12] with the conclusion that 


a 


-l 


1900 ± 50 ± 100 MeV for am q = 0.010 


( 6 ) 


( 1660 ± 110 ± 100 MeV for am, = 0.025. 

The first error is statistical and the second is our estimate of the systematic uncertainties 
implicit in choosing the S-P mass splitting as the physical quantity to set the scale. Thus 
we find for the mass of the a\ 
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(4) 


Z„ = (0|.4,|ai) = 

where f^ at is dimensionless. 

J 

The simulations were carried out on the Connection Machine CM-2 at the Supercom¬ 
puter Computations Research Institute at Florida State University. We used the ensemble 
of configurations generated by the HEMCGC collaboration with two flavors of dynamical 
staggered quarks [4], The configurations were generated using the hybrid molecular dy¬ 
namics (HMD) algorithm [5]. The size of the lattices is 16 3 X 32, the lattice coupling is 
f3 = 6 /g 2 = 5.6, and the dynamical quark masses in lattice units a are am q = 0.010 and 
0.025. Periodic boundary conditions were used in all four directions of the lattice. The total 
simulation length was 2000 simulation time units (with the normalization of ref. [6]) at each 
quark mass value. We analyzed lattices spaced by 20 HMD time units, for a total of 100 
lattices at each mass value. 

The spectroscopy was computed with six values of the Wilson quark hopping parameter: 
k = 0.1600, 0.1585, 0.1565, 0.1525, 0.1410, and 0.1320. The first three values are rather 
light quarks (the pseudoscalar mass in lattice units ranges from about 0.25 to 0.45), and the 
other three values correspond to heavy quarks (pseudoscalar mass from 0.65 to 1.5). Our 
inversion technique is conjugate gradient with preconditioning via incomplete lower-upper 
(ILU) decomposition by checkerboards [7]. For more details about the dynamical staggered 
fermion simulations see ref. [4,8]. Since we use sources for the propagators which are extended 
in space, we fix gauge to lattice Coulomb gauge using an overrelaxation method [9]. 

A few words are in order about our implicit assumptions. The state created by the 
correlator at asymptotic Euclidean time t is the on-shell (physical) oq, and so the matrix 
element calculated is the on-shell decay constant. In tau decay the produced oq is usually 
virtual (not at the peak of its resonance) and therefore is not “on shell”. In practice the 
off-shell production of the oq is entangled with its decay into three pions, which has to 
be modeled somehow. To make contact with tau decay, we follow the phenomenological 
analysis of Isgur, Morningstar, and Reader (IMR) [10] which parameterizes tau decay into 


3 



First seen as a resonance in pion-proton scattering, the axial vector meson cq remained 
an elusive prey due to its large width and the presence of strong background signals [ 1 ], 
Since 1986 the properties of the oq have been measured more precisely through the decay of 
the tau lepton to three pions: r —> v T a\ —> v T pir —> v T tttttt [ 2 ]. 

The decay of the oq is parameterized by the W —> oq vertex which, in the conventions 
of Tsai [3], is defined in terms of a dimensionful decay constant f ai as 

T^(W^ ai ) = -*/ 0l (p 2 )^. (1) 

In this definition f ai has units of [mass] 2 . The mass and decay constant are both quantities 
accessible to measurement from lattice QCD. In this work we describe a lattice calculation 
of the mass and decay constant of the oq meson and compare our results to experiment. 

The quantities measured on the lattice which give the oq mass and decay constant are 
two point correlation functions 

C,,(i) = (0|0,(i)0,(0)|0), (2) 

which, for large t, reduce to a single decaying exponential (plus boundary terms): 

ZjZj -aMt 

2aM 

where M is the mass of the lightest particle that couples to (9, a is the lattice spacing, and 
the Z’s are the appropriate source and sink matrix elements, Z{ = (0\Oi\M). In order to 
couple to the oq, the operator we use is the local axial-vector current A; = ^ 7 ^ 75 ^. For the 
operator which creates the state, we gauge fix to Coulomb gauge and use an axial current 
operator in which the quark and antiquark are created with Gaussian spatial distributions 
about a source point. To measure the decay constant we compute two kinds of correlation 
functions, ones in which the second operator is identical to the source and ones in which 
the second operator is the local current. A simultaneous fit to the two correlators allows 
an extraction of the mass and decay constant. We parameterize the matrix element for the 
local sinks as 
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Abstract 


We determine the mass and decay constant of the a\ meson using Monte 
Carlo simulation of lattice QCD. We find M ai = 1250 ± 80 MeV and f ai = 
0.30 ± 0.03 (GeV) 2 , in good agreement with experiment. 
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